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FORMATION OF LATE-TYPE SPIRAL GALAXIES: 
GAS RETURN FROM STELLAR POPULATIONS REGULATES DISK DESTRUCTION AND BULGE GROWTH 
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ABSTRACT 

Spiral galaxies have most of their stellar mass in a large rotating disk, and only a modest fraction in 
a central spheroidal bulge. This challenges present models of galaxy formation: galaxies form at the 
centre of dark matter halos through a combination of hierarchical merging and gas accretion along cold 
streams. Cosmological simulations thus predict galaxies to rapidly grow their bulge through mergers 
and instabilities, and to end-up with most of their mass in the bulge and an angular momentum much 
below the observed level, except in dwarf galaxies. We propose that the continuous return of gas 
by stellar populations over cosmic times could help to solve this issue. A population of stars formed 
at a given instant typically returns half of its initial mass in the form of gas over 10 billion years, 
and the process is not dominated by supernovae explosions but by the long-term mass-loss from low- 
and intermediate-mass stars. Using simulations of galaxy formation, we show that this gas recycling 
can strongly affect the structural evolution of massive galaxies, potentially solving the bulge fraction 
issue, as the bulge-to-disk ratio of a massive galaxy can be divided by a factor of 3. The continuous 
recycling of baryons through star formation and stellar mass loss helps the growth of disks and their 
survival to interactions and mergers. Instead of forming only early- type, spheroid-dominated galaxies 
(SO and ellipticals), the standard cosmological model can successfully account for massive late- type, 
disk-dominated spiral galaxies (Sb-Sc). 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: bulges 



1. INTRODUCTION 

The formation of disk-dominated gala xies r e mains a 
challenge for modern cosmology (White 2009; iBurkertl 
2009). Bright spiral galaxies with stellar masses around 
10 11 M at z = typical ly have bulge-to-disk mass ra- 
tios (B/D) of 0.2 to 0.4 (We inzirl et aLll2QQ9[ ). For in- 
stance, the Milky W ay itself has a B/D ratio of 0.35-0.40 
(jRobin et al.ll2003l ). Luminosity rati os in near-infrared 
and optical bands are even lower (jGraham fc Worlevl 
2008). However, cosmological simulations predict much 
higher bulge fractions, and consequently a baryonic an- 
gular momentum much lower than observed. The hi- 
erarchical assembly of dark matter halos drives suc- 
cessive galaxy mergers, followed by a rapid growth 
of central bulges. Indeed, a single major merger, 
even starting with a high gas fraction, will generally 
end- up with more baryons in t he bulge than in the 
disk (Springel & Her nquistl [20051 : iRobertson et al.l 120061 : 
Hop kins et al . 2009b). Successive minor mergers also 
drive bulge growth (Bournaud, Jog & Combes 2007). 

It has been realized recently that the assembly of 
baryons onto galaxies is far from being only driven by 
mergers. Rapid accretion of cold gas is a nother major 
mode of galaxy assembly (De kel et al.ll2009ah . This cold 
accretion mode could apparently feed disk-dominated 
galaxies, but it makes the disks so massive and turbulent 
that they violently fragment into giant clumps (Bour- 
naud & Elmegreen 2009; Dekel, Sari & Ceverino 2009; 
Agertz, Teyssier & Moore 2009; Burkert et al 2009). 
Clum p coalescence will rapidly fuel the bulge (jNoguchil 
1999; Bournaud, Elmegreen & Elmegreen 2007), and 
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signatures of this additional bulge formation process 
have been observed (jGenzel et al.ll2008t lElmegreen et al.l 
2009). Disk internal secular evo lution can also lead to the 
slow buildup of pseudo-bulges (jKormendv fc Kennicutt I 
l200l . 

A general result is that the rapid formation of mas- 
sive bulges appears unavoidable because of the com- 
bination of mergers and disk instabilities. The most 
disky galaxies in cosmological models have B/D around 
1 in the best cases and generally higher, in both the 
merger-driven and stream- fed dominant modes (e.g. , 
Governa toeTaIll2009al: IWhitd 120091 : iGibson et all 120091 : 
iScannapieco et aLl 120091 : Ceverino, Dekel & Bournaud 
2009). The bulge fraction problem is related to the 
well-known spin crisis: galaxies could gain angular mo- 
mentum from gravitational torquing by satellites (Por- 
ciani, Dekel & Hoffman 2002), but the early build-up 
of bulges quenches this process. Energy feedback from 
supernovae explosions can regulate star formation and 
keep the gas fraction higher, so that disk components 
survive mergers more easily. Stream-fed disks with ef- 
ficient supernovae feedback can then have realistic ro- 
tation velocitie s, in agreement with the observed Tully- 
Fisher relation (Mailer & Dekel 2002; Scannapieco et al. 
2008; Ceverino & Klvpin 2009; Governato et al. 2009a; 
IGibson et al. 20091: iPiontek fc Steinme tz 2009a), but the 
mass fraction in these fastly rotating disks remains too 
low compared to their bulges, meaning that the real 
angular momentum is low. High bulge fractions thus 
remain as a separate issue that supernovae feedback 
cannot solve, even when ma ny feedback models are 
tested in numerical simulatio ns ((Scannapieco et all 2009: 
IPiont ek & Steinmetz 2009b|). Apart for dwarf galaxies 
(jGovernato et al.l 2009b), the standard models appar- 
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ently fail to form high-spin, disk-dominated galaxies. 

A stellar evolution process that is more general than 
supernovae feedback is the continuous return of gas from 
stars of any mass, in particular through stellar winds 
and planetary nebulae (Faber & Gallagher! 1 1976f ). The 
gas return fraction of a coeval stellar population over the 
Hubble time can be as high as 50% for a Chabrier IMF, 
or around 40% for a Scalo IMF, with small dependencies 
on the metallicity (Jungwiert, Combes & Palous 2001; 
Lia, Portinari & Carraro 2002; Pozzetti et al. 2007). 
Supernovae account only for about 10% of this return 
fraction, which is dominated by low- and intermediate- 
mass stars. Half of the stellar mass forming at redshift 
z ~ 2 will thus be returned in the form of gas by red- 
shift 0, which can cool down in the disk and form new 
stars that will in turn return gas, in a continuous recy- 
cling process. This is another type of feedback, which 
is more a mass feedback than the energy feedback from 
supernovae explosions usually considered. This contin- 
uous mass return has alread y been included in a few 
cosmological simulat ions (e.g., iSemelin fc Com bes 2005; 
W iersma et al. 2009), but these simulations were limited 
to high redshift, or lacked resolution to clearly resolve 
disks and bulges, or the specific effects of gas return were 
not explored. 

Here we show that this known process has unsuspected 
effects on regulating the growth of bulges and the survival 
of massive disks in spiral galaxies. Using a cosmological 
simulation of the formation of a massive galaxy, which in 
the standard model ends- up with a B/D ratio larger than 
unity, we show that continuous gas recycling reduces the 
final B/D ratio to 40%, and the stellar light fraction in 
the bulge by a factor of 3. The mass and disk rotation 
speed of this galaxy are consistent with observed scaling 
relations, and its remarkably flat rotation curve is typical 
for late- type, high-spin spirals. The mass return of stars 
regulates the bulge growth and angular momentum dissi- 
pation in such proportion that a late-type spiral (Sb-Sc) 
is formed instead of an early-type galaxy. This shows 
how strongly stellar evolution can affect the formation of 
galaxies, and help to solve the issue of massive late- type 
spiral galaxy formation. 

2. SIMULATION 

Our model is based on a high resolution re-simulation 
of a galaxy identified in a large-scale cosmological sim- 
ulation. The initial cosmological run was performed as- 
suming standard A-C DM cosmolog ical parameters using 
the RAMSES code (jTevssierl [20021 ) and is described in 
iMartig et al.l (j2009l ). The chosen dark matter halo has a 
virial mass of 1.3xl0 12 M at redshift 0, in a field en- 
vironment. We avoid dense environments and halos un- 
dergoing many violent mergers so that the chosen halo is 
prone to hosting a disk-dominated spiral galaxy at z = 0. 
The mass assembly history of the re-simulated galaxy is 
shown on Figure [TJ and results from a combination of 
galaxy mergers and diffuse infall from cosmic gas fila- 
ments. 

Usin g the zoom-in technique presented in Ma rtig et al.l 
(|2009h . we re-simulate the evolution of the chosen galaxy 
from z = 5 to z = inside a 400 kpc-large zoom area. 
The spatial resolution (gravitational softening) is 150 pc 
and the mass resolution (particle mass) is 1.5 x 10 4 M 
for gas, 7.5 x 10 4 M for stars, and 3 x 10 5 M for 
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Fig. 1. — Stellar mass evolution of the studied galaxy as a func- 
tion of redshift. An intense phase of mergers takes place between 
z = 5 and z = 3. The subsequent evolution is much smoother, 
dominated by diffuse gas infall: the most important merger of this 
second phase takes place at z = 0.8 with a mass ratio of 6:1. 



dark matter. Gas dynamics is modeled with a sticky- 
particle scheme and star formation is computed with a 
Schmidt-Kennicutt law with an exponent of 1.5. The 
threshold for star formation is set to 0.03 M pc -3 , which 
corresponds to t he minimal densit y for diffuse atomic 
clouds formation ([Efmegreen 1 120021 ) . 

The galaxy is initialized at z = 5 as a very gas-rich 
disk, with a gas fraction of 0.5 with respect to total bary- 
onic mass. The initial dark halo follows a Burkert profile 
with a core radius of 5.9 kpc and a truncation radius of 
13.8 kpc. Gas and stars are in a disk with exponential 
scale-lengths of 540 and 230 pc, respectively. Starting the 
simulation at z = 5 captures 97% of the mass assembly 
history of this galaxy, as the total mass (including dark 
matter) at redshift 5 is only 3.7 x 10 10 M . For this rea- 
son and because numerous mergers rapidly destroy the 
proto-disk into a spheroid around which the regular bulge 
and disk components will gradually grow (Figure [2j) , our 
choice of starting with a disk is not influential in the final 
z = properties of the galaxy. 

We do not include the energy feedback from super- 
novae. There is no unique model for this, and feed- 
back is anyway most efficient in dwarf galaxies where th e 
structure is strongly influenced (iGovernato et a l. 2009b); 
the effects in mass ive galaxies are more moderate (e.g., 
lAgertz et al.ll2010l ). Not including feedback would be a 
serious problem at ver y high redshift (z > 5) becau se of 
the overcooling issue (jSommer-Larsen et al. I fl999). or 
in very high resolution (pc-scale) simulations to avoid 
the collapse of all gas into very dense clumps that would 
never be disrupted. This is not resolved in our simula- 
tions at resolutions of 100 pc or more. At these scales, 
feedback mostly regulates star formation, and we cali- 
brate the star formation efficiency in our simulations so 
that the gas consumption time matches typical observed 
values, so this effect of feedback is indirectly modeled. 

In turn, we implement the continuous gas return from 
high-, intermediate- and low-ma ss stars. To this a im, w e 
follow the scheme proposed by Jung wiert et al.l ((2001): 
each stellar particle represents a population of stars born 
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Fig. 2. — Time evolution of the studied galaxy. Left: large scale gas distribution at z = 5. The main galaxy is seen at the center; the 
box side is 400 kpc. Middle: evolution from z = 4 to z = of the face-on i-band surface brightness, with stellar mass loss (snapshot size 
50x50 kpc, log scale). Right: Edge-on views of the i-band surface brightness at z = in log scale. Each snapshot is 15x15 kpc, with 
contours at 17.25, 18.5, 19.75 and 21 mag arcsec -2 . In both simulations, a spheroid is first formed by successive mergers at z > 2. A 
large spiral disk is built at lower redshift, both from gas found in-situ and accreted from cosmic filaments and companion galaxies. Stellar 
mass loss supports disk survival in mergers and disk rebuilding after mergers, leading to a much lower B/D. At redshift 0, a late-type, 
disk-dominated Sb-Sc galaxy with a modest bulge is formed, instead of an early-type SOa-like system with a large and massive bulge and 
a low-mass disk. 



at the same time, whose return rate is given by: 

M = -M-^— (1) 
t-T v J 

where t is the age of the particle since its formation, 
T is set to 4.86 Myr to fit standard IMF return function 
(jJungwiert et al.f 2001). and C is such that the integrated 
return fraction of a stellar population over 10 Gyr is 40%. 
The returned gas is directly set cold, immediately avail- 
able for star formation. Real stellar winds would provide 
warmer gas, but overestimating the delay before it can 
form new stars would artificially increase the gas frac- 
tion and the capacity of disks to survive mergers: to 
avoid such favorable biases, we let the returned gas im- 
mediately available for star formation. 

We performed a "fiducial simulation" without gas re- 
turn from evolved stars, and a simulation with stellar 
mass loss. The structural properties were analyzed z — 
in both cases. We measured bulge-to-disk mass ratios 
with different techniques. First, we used azimuthally- 
averaged stellar density profiles: an exponential disk is 
fitted and the mass excess in the central regions corre- 
sponds to the bulge. We also performed a decomposition 
based on kinematics, and more precisely on the angular 
momentum of each star. We computed the total angular 
momentum of the gas disk in its inner 10 kpc, and set 
this as the z axis. For each stellar particle, we then com- 
puted e = jz/jcirc, the ratio of its angular momentum 
along the z axis to the angular momentum it would have 
if it were on a circular orbit (that is jd rc = r x v C i rc ). 
The distribution of the values of e shows two peaks: one 
around e = 0, that corresponds to bulge stars, and one 
around e = 1, that corresponds to disk stars. The limit 
between bulge and disk is set at the value of e corre- 
sponding to the minimum of the distribution between 
these two peaks. 

Finally, we also computed i-band surface luminosi- 
ties using the spectral evolutio n model PEGASE.2 
(|Fioc fc Rocca-Volmer angel Q999) to derive bulge-to- 



disk ratios in the z-band, assuming a solar metallicity: 
the i-band mass-to-luminosity does not vary much across 
the radial metallicity gradient of typical z = spiral 
galaxies. 

3. RESULTS 

The time evolution of the studied galaxy is shown 
on Figure [2j Interactions and mergers with companion 
galaxies rapidly transform the initial disk into a spheroid. 
At later times, star formation in gas left over after these 
early mergers and gas gradually brought in by cosmic 
gas flows and merging satellites builds a large disk, while 
interactions continue to grow a bulge - at a rate that de- 
pends on the model. At z = 0, both cases have a central 
bulge and an exponential disk with grand design spiral 
arms. A major difference, though, lies in the distribution 
of the z = stellar mass in the disk and bulge compo- 
nents. 

Without stellar mass-loss, the bulge to disk mass ra- 
tio derived from the stellar mass profile is slightly above 
1, which is typical for massive "disk" galaxies in current 
A-CDM models (see Introduction). The disk component 
contains barely half of the stellar mass, even though we 
include both the thin and thick disks. When the con- 
tinuous stellar mass loss is included, the bulge to disk 
mass ratio is reduced to 0.40: the bulge is less massive, 
with a lower density and a smaller size (right panels of 
Figure [2j) , while the disk is denser with a similar scale 
length (Figure [3] and Table QJ . 

Other definitions of B/D lead to the same conclusion. 
A kinematical decomposition even shows a more drastic 
effect, with B/D decreasing from 1.40 to 0.32 (Figure |4] 
and Tabled]). In the i band luminosity, B/D decreases 
from 0.49 to 0.16. Similar values were obtained when 
the simulation was run with a twice lower spatial reso- 
lution and a eight times higher particle mass, suggesting 
reasonable convergence of the simulation. 

The effect of stellar mass loss on the z = struc- 
tural properties of the galaxy is thus major, although 
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TABLE 1 

Disk and bulge properties at z = with and without stellar mass loss 



Fiducial run With mass loss 



Stellar disk mass (Mq) 


6.1 x 10 10 


8.8 x 10 10 


B/D (stellar mass profile) 


1.01 


0.40 


B/D (kinematics) 


1.40 


0.32 


B/D (i band) 


0.49 


0.16 


Disk exponential scale length (kpc) 


4.5 


5.0 


Optical radius R25 (kpc) 


22.0 


19.8 


Bulge half-mass radius (kpc) 


1.3 


1.6 


Absolute / band magnitude 


-22.8 


-22.9 


HI linewidth (km s _1 ) 


583 


591 



10 4 




5 10 15 20 25 

Radius [kpc] 



Fig. 3. — Stellar surface density profiles for the fiducial simulation 
and the one including stellar mass loss. Both galaxies have a central 
bulge and an exponential disk, but their relative importance differs. 
On the inset is shown the ratio between the stellar surface density 
when stellar mass loss is included to the surface density in the 
fiducial simulation: stellar mass loss reduces the bulge density and 
increases the disk one, thus decreasing B/D. Differences are seen 
even on log- scale radial density profiles. 

we studied quite a massive galaxy for which preserving 
disks is harder than in lower-mass systems. Continuous 
stellar mass loss transforms a bulge-dominated system, 
typically classified as SQ/ Sa in the Hubble classification 
([Graham fc Worlevll2QQ8h into a disk-dominated Sb/Sc- 
like galaxy, with a bulge fraction typical for the Milky 
Way and the majority of massive spirals in the nearby 
Universe. 

As a result from the different mass distribution, the 
rotation curve at z = has a different profile (Figure |4]). 
In the fiducial case, the rotation curve is peaked a low 
radii before decreasing and reaching a flat plateau. This 
is typical for bulge-dominated, SO-Sa galaxies (e.g., M31 
- Chemin Carignan fc Foste r 2009; NGC 3031, 4736, 
or 2841 in lde Blok eraT1l2009f ). The regulation of bulge 
growth by continuous stellar mass-loss results in a flater 
rotation curve, which is typical of disk- dominated spira l 
galaxies like NGC 2903, 3198, 3621 (jde Blok et al.ll2009l ). 
The comparison of rotation curve with typical cases also 
corresponds to a shift by two categories in the Hubble 
classification scheme, from early- type disky galaxies to 
late- type disk-dominated galaxies, in agreement with the 
direct estimates of bulge fractions. 

We also derived the equivalent of an observed atomic 
gas (HI) linewidth: we select gas in regions denser than 
one atom per cubic centimeter in order to include the 
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Fig. 4. — Kinematics in the two simulations. The top panel 
shows the rotation curves. The prominent bulge in the fiducial 
run causes a central bump in the circular velocity profile, while 
this profile is much flatter for the simulation including stellar mass 
loss. On the inset is shown the radial velocity distribution for the 
galaxy seen edge-on in both simulations: the full width at half 
maximum is used for comparisons with observations of the Tully- 
Fisher relation. The bottom panel shows the distribution of the 
ratio of the angular momentum of the stars along the z axis to the 
angular momentum they would have on circular orbits: bulge stars 
are found around e = while disk stars are found around e = 1. 
The bulge component is significantly reduced under the effect of 
stellar mass loss. 

warm/cold gas disk without the hot halo, simulate the ve- 
locity spectrum for the galaxy, and compute its FWHM 
averaged on several edge-on projections. The / band 
absolute magnitude to HI linewidth ratio (see Table [1]) 
puts both systems very close to the observed Tully- 
Fisher relation for local spiral galaxies as observed by 
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ISpringob et al. I ([2007D . The Tully-Fisher relation is nev- 
ertheless not the best direct estimate of the baryonic an- 
gular momentum of a galaxy since it measures the disk 
rotation speed at large radii without indicating whether 
or not a large fraction of the stellar mass lies in a con- 
centrated, low angular momentum bulge. 

A similar effect has been also found in two other sim- 
ulations. We measured bulge and disk mass follow- 
ing the decomposition based on kinematics. In one 
case (a galaxy with M* = 8 x 10 9 M at z = 5 and 
1.3xlO n M at z = 0), B/D decreases from 1.24 to 
0.75 when mass loss is included, in the other case (with 
M* = 5.3 x 10 9 M at z = 5 and 4 x 10 10 M at z = 0), 
it decreases from 0.94 to 0.49. 

4. CONCLUSION 

Our simulation of cosmological galaxy formation shows 
that the continuous gas return from high-, intermediate- 
and low-mass stars over cosmic times can strongly reg- 
ulate the growth of bulge and the dissipation of angular 
momentum in disk galaxies. There are at least two ways 
in which stellar mass loss can support the survival of 
large, massive disks. First, it can keep the gas frac- 
tion higher in young galaxies before a merger occurs, 
which helps disks surv ival against destruction into bulges 
(Hop kins et al.ll20 09a). Second, the bulge, stellar halo 
and thickened disk components release fresh gas after 
interactions and mergers: this decreases their mass and 
increases that of the thin gas disk. The disk re-formed by 
gas returned by a bulge or halo should have a low angular 
momentum, but can be torqued by satellites into a large, 
high-spin disk that gradually forms new stars (an exam- 
ple of such an interaction is shown at z = 1 on Figure [2]) . 
The radial migration of disk stars ca n also arise from res - 
onant interactions with spiral arms: Roskar et al.l {2008) 
have shown, in a Milky- Way like galaxy, that up to 50% 
of stars in the solar neighborhood could have migrated 
from the inner disk. In our model, the large disk scale 
length is preserved by the combination of stellar mass- 
loss and radial migration following external interactions 
and/or internal evolution. 

Overall, the effect of stellar mass loss on disk survival 
and regulation of bulge growth is major. It reduces the 
bulge-to-disk ratio by a factor up to ~ 3 in both the 
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stellar mass and light. This factor is what was typically 
missing in cosmological models to account for the forma- 
tion of massive, disk-dominated, late-type spiral galaxies. 

This does not imply that all galaxies would end-up 
with a spiral-like morphology at redshift zero, and we 
briefly presented cases of galaxies with other merger 
and star formation histories where the decrease of B/D 
is weaker. Gas return promotes disk survival against 
mergers and disk instabilities, but will preserve a disk- 
dominated galaxy only if the mergers are not too numer- 
ous, do not happen too late, and the internal instabili- 
ties not too violent. For instance, a galaxy undergoing 
a major merger at low redshift would still end-up as a 
z = elliptical or lenticular, because of a stellar pop- 
ulation globally too old to return a significant amount 
of gas. Such z = ellipticals are indeed found in our 
co smological simulati ons with the same technique (e.g. 
in lMartig et al.l 12009). Gas consumption, stripping and 
strangulation in groups, will also support the survival 
of early-type elliptical and SOs. Thus, gas return from 
stellar populations can explain the origin of late-type 
galaxies, without challenging the formation of early- type 
galaxies in systems that undergo later mergers or more 
violent internal instabilities. 

Other sorts of feedback processes, like the energy re- 
leased by supernovae or the radiation pressure from 
young massive stars (Murray, Quataert & Thompson 
2009) could further regulate the bulge growth especially 
in low-mass galaxies, potentially forming almost bulge- 
less galaxies. However, we propose that the formation of 
massive late- type galaxies in the A-CDM Universe is ex- 
plained not by an unknown combination of minor factors, 
but mostly by one major effect, namely by the gas return 
from stellar populations all across their initial mass func- 
tion. This process can largely help to solve one of the 
main challenges in the origin of modern galaxies, namely 
the origin of late-type, disk-dominated spiral galaxies like 
the Milky Way. 



We acknowledge support from GENCI-CCRT (Grant 
2009-042192) and Agence Nationale de la Recherche 
(ANR-08-BLAN-0274-01), and help from Romain 
Teyssier with the re-simulation technique. 



Dekel, A., Sari, R., & Ceverino, D. 2009b, ApJ, 703, 785 
Elmegreen, B. G 2002, ApJ, 577, 206 

Elmegreen, B. G., Elmegreen, D. M., Fernandez, M. X., & 

Lemonias, J. J. 2009, ApJ, 692, 12 
Faber, S. M., & Gallagher, J. S. 1976, ApJ, 204, 365 
Fioc, M., & Rocca-Volmerange, B. 1999, |astro-ph/9912179l 
Genzel, R., et al. 2008, ApJ, 687, 59 
Gibson, B. K. et al. 2009, IAU Symposium, 254, 445 
Governato, F., et al. 2009a, MNRAS, 398, 312 
Governato, F., et al. 2009b, Nature, 463, 203 
Graham, A. W, & Worley, C. C. 2008, MNRAS, 388, 1708 
Hopkins, P. F., Cox, T. J., Younger, J. D., & Hernquist, L. 2009a, 

ApJ, 691, 1168 
Hopkins, P. F., et al. 2009b, MNRAS, 397, 802 
Jungwiert, B., Combes, F., & Palous, J. 2001, A&A, 376, 85 
Kormendy, J., & Kennicutt, R. C, Jr. 2004, ARA&A, 42, 603 
Lia, C, Portinari, L., & Carraro, G. 2002, MNRAS, 330, 821 
Mailer, A. H., & Dekel, A. 2002, MNRAS, 335, 487 
Martig, M., Bournaud, F., Teyssier, R., & Dekel, A. 2009, ApJ, 

707, 250 



6 



Martig & Bournaud 



Murray, N., Quataert, E., & Thompson, T. A. 2009, ApJ, 709, 191 
Noguchi, M. 1999, ApJ, 514, 77 

P iontek, F., fe Stei nmetz, M. 2009a, MNRAS submitted, 

arXiv:0909.4156 
P iontek, F., fe Stei nmetz, M. 2009b, MNRAS submitted, 

arXiv:0909.4167 
Porciani, C, Dekel, A., & Hoffman, Y. 2002, MNRAS, 332, 325 
Pozzetti, L., et al. 2007, A&A, 474, 443 
Robertson, B. et al. 2006, ApJ, 645, 986 

Robin, A. C, Reyle, C, Derriere, S., & Picaud, S. 2003, A&A, 
409, 523 

Roskar, R., Debattista, V. P., Quinn, T. R., Stinson, G. S., & 
Wadsley, J. 2008, ApJ, 684, L79 



Scannapieco, C, White, S. D. M., Springel, V., & Tissera, P. B. 

2009, MNRAS, 396, 696 
Scannapieco, C, Tissera, P. B., White, S. D. M., & Springel, V. 

2008, MNRAS, 389, 1137 

Semelin, B., & Combes, F. 2005, A&A, 388, 826 
Sommer-Larsen, J., Gelato, S., & Vedel, H. 1999, ApJ, 519, 501 
Springel, V., & Hernquist, L. 2005, ApJ, 622, L9 
Springob, C. M., Masters, K. L., Haynes, M. P., Giovanelli, R., & 

Marinoni, C. 2007, ApJS, 172, 599 
Teyssier, R. 2002, A&A, 385, 337 

Weinzirl, T., Jogee, S., Khochfar, S., Burkert, A., Sz Kormendy, J. 

2009, ApJ, 696, 411 

Wiersma, R. P. C, Schaye, J., Theuns, T., Dalla Vecchia, C, &; 

Tornatore, L. 2009, MNRAS, 399, 574 
White, S. D. M. 2009, IAU Symposium, 254, 19 



